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Abstract 

By adjusting the parameters of the new phenomenological nuclear potential introduced recently 
by P. Salamon and T. Vertse [ll (SV) to the global potential sets of Percy and of Becchetti and 
Greenlees, it is shown that the SV potential well approximates the cut-off Woods-Saxon (CWS) 
form even if its second term is set to zero. The important difference between the CWS and SV 
potentials is that CWS jumps to zero at the cut-off radius, while the SV potential goes to zero 
smoothly. We fitted the single-term SV potential to the global potentials, and compared the 
spectra and the pole trajectories of neutron single-particle states generated by some light neutron- 
deficient nuclei that are most often used in positron emission tomography (PET): ^^N, ^^O and 
^*F. The CWS and the SV energies agree reasonably well, and only the pole trajectories, which 
span extremities, are substantially different. 
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The positron emitting nuclides used in PET are usually produced by irradiation of light 
targets by light charged particles accelerated by cyclotrons. The most productive reactions 
are -'^'^N(p, a)-'^^C, ■^^C{p,ny^l<l, ^^N(p, n)-'^^0, and ^^0(p, n)^^F. In three of the four reactions 
listed above neutrons are emitted, so the exit channel is neutron scattering from ^^N, ^^O, and 
-^^F, respectively. Neutron scattering is generally described in the framework of the optical 
model, in which the mean field of the target nucleus is represented by a phenomenological 
optical potential. 

The geometry of the optical potential terms is generally fixed and only the strengths of 
its terms are assumed to be energy-dependent. The optical potential is constructed almost 
exclusively from cut-off Woods-Saxon (CWS) terms and their derivatives. The strict finite 
range (SFR) is emphasized j^] since whenever the radial Schrodinger equation is solved 
numerically, we in fact use SFR potentials since we match the numerical solution to a linear 
combination of the incoming and outgoing Coulomb waves at a finite distance r = -Rmatch, 
and calculate the S'-matrix from this matching. For neutron scattering the Coulomb waves 
reduce to the solutions of the Riccati-Hankel differential equation. 

The CWS can be written as: 

\/^wS(r, R, a, i?„,ax) = -K)/cws(r, R, a, R^,^) , (1) 
where the radial shape is 

, if r < i?max , , 

/cws(r,/2,a,/2max) = < (2) 

I , if r > i?max • 

Nuclear scattering may be described in terms of resonances. It was shown earlier [L] that 
in the CWS potential the positions of broad resonances do depend on the value of the cut-off 
radius i?max [21, |3[, although it has no physical meaning. This is caused most probably by 
the discontinuity of the CWS form at Rmax- Therefore, the cut-off radius is an important 
parameter of the CWS form in Eq. ([2]) . The two-other parameters of the CWS form are the 
radius R and the diffuseness a. 

In Ref. [l| two of us proposed a more convenient SFR potential, the Salamon-Vertse (SV) 
potential instead of the CWS potential. The SV potential tends to zero smoothly, and even 
its derivatives become zero at and beyond a finite radial distance. Thus the SV potential 
has the attractive mathematical property that its derivative exists in the whole r G (0, oo) 
range. 



The SV potential may be specified [ij as a linear combination of the function 

/(r,p)=r"^ * " (3) 

, if r > p , 



and its derivative with respect to r, /'(r, p). The SV potential may thus be written [3| as 
a product of a potential depth Vq and a radial shape fsyi^r, c, po, pi) with three adjustable 
parameters po, pi and c: 

\/SV(^) = -K)/sv(r,c,po,pi) , (4) 
fsY{r,c,po,pi) = f{r,po) - cf'{r,pi) . (5) 

To describe attraction, the coefficient c has to be non-negative. With the usual choice 
(po > Pi), the SV form vanishes at po, i.e., /sv(r > po,c, po,pi) = 0. The SV radial form 
has three parameters just like the radial form of the CWS potential. 

To make the SV potential resemble the CWS shape as much as possible, we fit its three 
parameters to the radial form factor of the CWS shape /ws by minimizing 

A(po,Pi,c)= / [/sv('",c,po,pi) - /ws('^,^,a,^max)] dr. (6) 
Jo 

The integration in Eq. can be performed by using equidistant mesh-points = ih over 
the range of the integration: 



A(po,Pi,c) = ^ [/sv('^i,c,po,pi) - /ws('"i,^,a,^max)]^ • (7) 



[/svl'^i,C,Po,Pij - /wsl'"i,^,a,^maxj '^ 
i=l 

There are several families of global optical model potentials with CWS real parts. In 
most cases the radial shapes (the radius parameter R and diffuseness a) are fixed, and only 
the strengths Vq are varied. The spin-orbit potential for a particle with a spin of s = is 

V;^^(r, i?so, aso, i?max) = K^^Wslr, Rso, aso, ^max) 2{t- s) , (8) 

with a radial form involving the derivative of the central potential 

r 

The radial factor of the spin-orbit term is also cut to zero at -Rmax (the derivative /cws('") 
would be infinite at i?max)- It is easy to define the spin-orbit potential for the SV potential 
since f'svij") exists everywhere: 

(r, c, Po, pi) = V!^hsv{r, c, po, pi) 2{l- s) , (10) 



with 

hsv{r, c, po, pi) = --f'svir, c, po, pi) . (11) 



r 



1 /3 

The radii of the potential terms for a target with mass number are scaled by Ar^. , i.e., 

1 /3 

R = tqAc^ . The radii and diffuseness can be different for the different terms of the potential. 



There are two cWca, global nucleon potential set. given by Perey and by Becchetti 

and Greenlees [5| long ago. In a more recent study of the isospin-dependent global neutron- 



nucleus optical model potential jo], the old Becchetti-Greenlees potential was used as a 
starting point. In this work, we also use the Perey and the Becchetti-Greenlees parameters 
for simplicity. Since the SV potential is quite new, and its mass number dependence has not 
been explored [l|, we fitted the SV parameters to the global potential shapes by minimizing 
Eq. ([6]) for a broad range of mass numbers. The A^ dependence of the resulting SV parame- 
ters is quite regular except in the region of the light nuclei, and the SV potentials reproduce 
the CWS shape quite well. The mixing coefficient c decreases with decreasing At as one can 
see in Fig. [H To avoid the discrepancy for light nuclei, we imposed the requirement that the 
derivative of the SV potential should also be similar to the derivative of the CWS shape by 
minimizing 

dfsY in , ci , Po , Pi ) dfws in) 



A(po, Pi, ci) = ^ [/sv(ri, ci, Po, Pi) - fwsiri)] + X 



i=l 

.2 



(12) 



dr dr 

With a value of A = 25 fm^, the resulting SV potential became reasonably smooth and 
similar to the CWS shape that we wanted to approximate. 

With the Perey j4| parameters for the central potential (ro = 1.25 fm and a = 0.65 fm) 
we obtain that the best-fit SV parameters are po = l.'ShA!^'^ fm, pq — pi = 3.2a fm (a the 
diffuseness), c = —0.03 + 0.004At. With the Becchetti-Greenlees (5[ geometry (ro = 1.17 
fm and a = 0.75 fm), the radius of the central potential is, of course, somewhat smaller, 
but the best-fit SV parameters behave similarly to the case of Perey's potential. With Aj' 
decreased, the mixing parameter c decreased for both global potentials as one can see in 
Fig. [U The second term is apparently needed to produce the nearly constant inner section 
of the potential, and since the inner part of a light nucleus is very small, the second term has 
no significance for them. For example, the SV parameters conforming the Perey geometry 
are po = 5.084 fm, pi = 3.244 fm and c = 0.040, while those reproducing the Becchetti- 
Greenlees shape are po = 4.957 fm, pi = 2.728 fm and c = 0.011. Thus, for the light nuclei 
under consideration, it is reasonable to take simply c = 0. 



TABLE I: Geometrical parameters of the SV and CWS potentials for N, O and F nuclei. All 
distances are in fm units. 

At po R a 

4.79 2.94 0.65 
1^0 5.01 3.08 0.65 
i^F 5.3 3.28 0.65 



In Fig. m the CWS shape and its derivative with the SV potential and its derivative 
(c = 0) are compared for the Perey geometry apphed to the ^^F target. The CWS parameters 
are Hsted in Table [B This shape is described best by a c = SV form with a range po = 5.3 
fm. The SV shape is similar to the CWS shape although now we have a single parameter 
to adjust. The best-fit value of p is shifted a bit from the optimal value of po found with 
two terms (from po = 4.975 fm to p = 5.3 fm). The difference between the CSW and the 
SV potentials grows with the distance, as the SV shape goes to zero much faster than the 
CWS shape. 

It is interesting to see how these differences influence the single-particle energies in the 
two potentials. In Table Hi] we show the neutron single-particle energies e„/j calculated in 
the nucleus ^^F with CWS potential parameters = 60 MeV, tq = 1.25 fm, a = 0.65 
fm, = 7 MeV and with the fltted SV potential. The spin-orbit potential terms were 
calculated according to Eqs. (|8]) and (fTOj) . respectively. The strength of the SV potential 
has been readjusted so as to reproduce the Is energy obtained with the CWS shape, which 
resulted in Vq = 59.934 MeV. In Fig. [2] one can observe that the difference between the 
original Perey form and the adjusted SV form is magnified in the derivatives, which are used 
in the spin-orbit potentials h-wsii^) and hsy{r). The spin-orbit strength for CWS potential 
was Vso = 7 MeV. This value caused a gap of 1.93 MeV between the Op orbits and a gap of 
3.2 MeV between the Od orbits. With a V^^^ = 30 MeV, we were able to produce roughly 
the same gaps (1.97 MeV for / = 1 and 3.28 MeV for / = 2). 

We obtained similar agreements for the two other nuclei, ^^O and ^^N. Table UTTl represents 
the calculated single-particle energies for ^^N. 

We can conclude that for light nuclei the c = version of the SV potential is a good 
approximant to the phenomenological CWS potential, since it produces almost the same 
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TABLE II: Neutron single-particle energies (in MeV) in the CWS and the SV potentials outside 
the ^^F core. 



n,l,j Energy (CWS) Energy (SV) 



OSl/2 


-38.926 


-38.119 


0p3/2 


-23.998 


-23.610 


Opi/2 


-22.0674 


-21.640 


lSl/2 


-7.697 


-7.697 


0^5/2 


-8.985 


-9.048 


04/2 


-5.779 


-5.769 



TABLE III: Neutron single-particle energies (in MeV) in the CWS and the SV potentials outside 
the ^^N core. 



n,l,j Energy (CWS) Energy (SV) 



Osi/2 


-35.045 


-36.746 


0p3/2 


-18.620 


-20.368 


Opi/2 


-16.318 


-17.958 


l'Sl/2 


-3.400 


-3.400 


04/2 


-3.067 


-4.247 


04/2 


- 0.003 


- 0.548 



spectrum as the CWS potential. 

To see whether the SV potential is a good substitute for the CWS form, one has to 
compare their performances for positive energies as well. It is, however, difficult to compare 
continuum entities. What we do instead is a comparison in terms of 5'-matrix poles, whereby 
we can also check whether the features of the CWS potential observed [2] for heavy nuclei 
are present for light nuclei as well. It has been found in Ref. {2] that the real parts of the 
starting points of the s-wave pole trajectories follow Newton's rule: i?Re(A;„) fti nvr 3]- Here 
R denotes the range of the potential beyond which it is zero. In the CWS potential, R is 
the cutoff radius i?max- The Re(A;„) values calculated for the CWS potential follow this rule 
approximately if n is not very small. 



6 



The high-precision numerical solution of the radial equation was performed by using 
Ixaru's method [B]- In Fig. [3] the pole trajectories are shown for ^^F in a CWS potential 
well with -Rmax = 15 fm. By fitting the calculated Re(A;„) values for n < 5 with a first-order 
polynomial, we can deduce R = 14.8 fm. This agrees reasonably well with the R^ax = 15 
fm value used. 

In Fig. m we show the calculated Re(A;„) values for the corresponding SV potential. If 
Newton's rule were valid, the trajectories would start from the vicinity of these lines. The 
starting points, however, overshoot the trajectories. Thus the SV trajectories satisfy the 
rule less accurately. From the starting points of n > 3 we obtain R = 4.84 fm for the slope. 
The deviation from the pq = 5.3 fm value is not negligible and even a bit larger than the 
deviation found in Ref. It is more important that the pole trajectories produced by the 
SV potential behave quite regularly in contrast to the CWS potential in which the n = 1 
trajectory shows a strange circulating path. Similar strange behavior was also found for 
^°®Pb {q] and was attributed to the sharp cut-off of the CWS potential. 

The substantial differences between the pole trajectories may cause some concern as to 
the performance of the SV potential in scattering calculations, but the trajectories span a 
broad range of potential depths, which extend to physically irrelevant regions. From the 
point of view of scattering calculations in the region of isotope production, the most relevant 
property of the potentials is their low-energy behavior near the threshold. By putting some 
of the states at the threshold, we could also show that depths of the Percy and the SV 
potentials are almost the same, which demonstrates that the two potentials must behave 
very similarly for low-energy scattering relevant to the production of PET isotopes. 

This work was supported by the ENIAC CSI Project (No. 120209). 



[1] P. Salamon and T. Vertse, Phys. Rev. C 77, 037302 (2008). 

[2] J. Darai, A. Racz, P. Salamon, and R. G. Lovas, Phys. Rev. C 86, 014314 (2012). 

[3] A. Racz, P. Salamon, and T. Vertse, Phys. Rev. C 84, 037602 (2011). 

[4] F. G. Perey, Phys. Rev. 131, 745 (1963). 

[5] F. D. Becchetti and G. W. Greenlees, Phys. Rev. 182, 745 (1969). 

[6] Xiao-Hua Li, Lie- Wen Chen, Nucl. Phys. A874, 62 (2012). 



7 



[7] R. G. Newton, Scattering Theory of Waves and Particles, Springer Verl., New York Inc., 1982. 
[8] L. Gr. Ixaru, Numerical Methods for Differential Equations and Applications, D. Reidel, Publ. 

Comp., Dordrecht /Boston/Lancaster, 1984. 
[9] T. Vertse, R. G. Lovas, P. Salamon, and A. Racz, talk on the Int. Conf. Nuclear Structure and 

Dynamics II, held in Opatija, Croatia, July 13-19, 2012, the Proceedings will be published by 

the AIP later. 



8 




FIG. 1: Dependence of the mixing coefficient c on tlie target mass- number ^t- 
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r [fm] 

FIG. 2: Radial shapes of the CWS and SV (c = 0) potentials and its derivatives for ^^F. 
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FIG. 3: Trajectories of the poles of the ^-matrix in Z = partial wave in the CWS potential of 
Perey's shape. The vertical lines show the Re(A;) = values. 
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FIG. 4: Trajectories of the poles of the S'-matrix in Z = partial wave in the SV potential 



Perey's shape. The vertical lines show the Re(A;) = ^ values. 
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